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Surface solitons in VT symmetric potentials 
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We investigate light beam propagation along the interface between linear and nonlinear media 
with parity-time (VT) symmetry, and derive an equation governing the beam propagation. A novel 
class of two-dimensional VT surface solitons are found analytically and numerically, and checked 
to be stable over a wide range of VT potential structures and parameters. These surface solitons 
do not require a power threshold. The transverse power flow across beam is examined within the 
solitons and found to be caused by the nontrivial phase structure of the solitons. VT surface solitons 
are possibly observed experimentally in photorefractive crystals. 
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Surface solitons were found to be able to propagate 
along the interface between two media with real index 
of refraction, and showed many interesting properties 
which had no analogues in single medium alone [l|43[- 
It was also reported that surface solitons could exist at 
the interface between two media where at least one of 
them is nonlinear. Theoretical and experimental stud- 
ies on the surface solitons have been reported in two- 
dimensional domain where many interesting phenomena 
associated with nonlinear surface waves were predicted 
and observed [J-lZl 

On the other hand, Bender and Boettcher found a non- 
Hermitian Hamiltonian might have an entirely real eigen- 
value spectrum [8[ when potential field met parity-time 
(VT) symmetry, which implied eigenstate might exist in 
VT symmetric complex potential field. In optics, com- 
plex potentials mean a gain and loss contained in optical 
media. Consequently, the solitons in VT symmetric op- 
tical potentials were explored extensively in theory and 
experiments thereafter. Mussilimani studied VT soli- 
tons in one-dimensional (ID) and two-dimensional (2D) 
media [9|, LL0[. The basic properties of Floquet-Bloch 
modes in VT symmetric optical lattices were examined 
in detail [ll|. In 2010, optical VT symmetry was ob- 
served experimentally in an optical coupled system in- 
volving a complex index potential [12[. VT symmetry 
breaking was also demonstrated in complex optical po- 
tentials 13, LL4J . The effect of localized modes in lattices 



of size N with exponentially fragile VT symmetry was 
also studied [15]. The light transmission in large-scale 
temporal VT symmetric lattices was studied [16[ . It was 
also found that gray solitons might exist and be stabilized 
in VT symmetric potentials [17[. However, all investiga- 
tions mentioned above dealt with the VT symmetry and 
VT solitons only within one complex index medium. The 
surface solitons at the interface between two complex in- 
dex media, a generalized surface solitons, have not been 
explored. It is still unknown if such surface solitons may 
exist and stabilize. 



In this letter, we explore such surface solitons in two 
complex index media composed of self-focusing and lin- 
ear media, and find that the 2D surface solitons may 
exist and stabilize at 2D interface when the optical po- 
tential of the complex index self-focusing medium meets 
VT symmetry and the interface is oriented in a particu- 
lar direction. We gain the 2D VT surface solitons, and 
test their stability numerically. 

We express the optical field as E(x,y,z,t) — 
E(x,y)e l ( Tz ~ UJt \ the refractive index in nonlinear and 
linear media as noi + n r i (x, y) + inn (%,y) + n 2 1 E | 2 and 
^02 + n r 2(x,y) + ini2(x,y) : respectively, where u is the 
frequency, T is the wave vector, noi and no2 are the back- 
ground refractive indices, n r i, r 2(x,y) is the real index 
profile and nn,i2(x,y) the imaginary part (gain or loss) 
oiVT symmetric potentials, and n2 is nonlinear Kerr co- 
efficient [9|, llO( • Substituting the field and the refractive 
index into the Helmholtz equation, we get, 



dz 



d 2 E d 2 E 



(k 2 n 2 01 - t 2 )e 



dx 2 dy 2 

+2fco[n i(n r i + inn) + n in 2 \E\ 2 }E 



0, 



(1) 



in the nonlinear medium when the conditions, fco = cj/c, 
n r \ i\ « noi and n2\E\ 2 << noi, were met, and 



dz 



d 2 E d 2 E 
dx 2 dy 2 

+2/CQn 2(^r2 



(kin, 



02 



-T 2 )E 
in i2 )E = 0, 



(2) 



in the linear medium as n r 2,i2 << ^02 • 

Introducing the scaled dimensionless quantities as 
C = z/(2Twl), £ = x/w , n = y/w , q = 
n§ 1 /[konoiWo(2n2)z]E, R = 2noiklwl(n r i + inn) and 
R = 2no2kQWQ(n r 2 + m^), respectively in nonlinear and 
linear media, where wq is an arbitrary scaling factor, 
Eqs.(pp) and ([2]) can be changed as, 



.dq 



Vlq + Rq + G iq +\q\ 2 q = 0, 



(3) 



i-Tf + Vl<? + Rq + G 2 q = 0, 



(4) 



where G x = w 2 [(k n 01 ) 2 - T 2 ], G 2 = w 2 [(k n 02 ) 2 - T 2 ] 
and Vi = d 2 /d£ 2 +d 2 /dr] 2 . Setting q = q'e iG ^, Eqs.© 
and dU can be changed as follows, 



.f + vl,< 



Rq' 



\q'\ 2 q> 



0, 



i^r + V 2 ± q' + Rq' + AGq' = 0, 



(5) 



(6) 



with the propagation constant 6 = 4. 

To shed light on the formation of the VT surface solu- 
tion, we examine the transverse power-flow density across 
the beam, which can be described by Poynting vector [9[ , 
S = (i/2)(uVu* - u*\7u), where V = d/d^ + id/drj, and 
intuitively indicates the energy exchange. For the soliton 
in Eq.(12), we get, 



where AG = G2 — G\ = ^0^0 ( n 02 — n oi)- Based on the 
step function U(x,y) which is equal to 1 or 0, respectively 
in nonlinear and linear media, Eqs.(j5]) and (J6j) can be 
integrated into one equation as follows, 



•¥ +vi «' 



Rtf + (l-U)AGq f + U\</\ 2 q' = 0, (7) 



Assuming a stationary soliton solution to be in th 
form q'i^rjX) — u {^-> r \) e%hCi (iSl? where u is a comple 
function and b is the propagation constant, Eq.© ca 
now be deduced to, 

V\u + [R + (1 - U)AG]u + U\u\ 2 u -bu = 0. (8 

R + (1 — U)AG is a VT symmetric potential becaus 
R has VT symmetry. As a result, Eq.® describes 
stationary VT soliton [9] in the VT potential, R + (1 - 
U)AG. 

To seek for possible analytical stationary soliton sc * 
lutions to Eq.(j8]), a special VT symmetric potentia 
R = V + iW ", is constructed, which is expressed by, 



w, 



V = (4 + ^) sec/i 2 (£ + r?) + 4sec/z 2 (£ 
- [1 + e a ^-^]- 2a sec fr 2 (£ + 77) sec fr 2 (£ ■ 

W = W sec ft(f + rj) tanh(£ + 77), 



-77) 



(9) 
(10) 



(1 - C/)AG = 



2a 4 +2a^ 



2a 3 -4a 2 tanh(£-77) 
l+e«(»7-€) 



,(n) 



where the step function E/(£, 77) is approximated by the 
function, l/[l + e Q; ^ _77 ^, to guarantee its differentiability. 
Meanwhile, this approximate step function also implies 
that the interface between nonlinear and linear media 
coincides with the line £ — 77 = 0. VFo is the amplitude of 
the imaginary part of the VT symmetric potentials, a is 
a positive factor which describes the degree of proximity 
between 1/[1 + e a ^~^} and U(£,rf). 

For the VT symmetric potential described by Eqs.(|9])- 
(pjj, it is possible to realize it in a photorefractive crystal 
by means of appropriate light pumping through specially 
designed masks, as shown in [12|, and an analytical sta- 
tionary soliton solution to Eq.flSJ) is found as follows, 



sec/i(£ + 77) sec/i(£ — 77) ^wj 



[l^e a ^~^} a 



—jP- arctan[sinh(£+?7)] 



,(12) 



X _ (I-H)Wq sec h 2 (j-ry) sec/i 3 (£+r?) 

D - 6 [l + e «(£-r ? )]2« 



(13) 



It is interesting to note that the energy-flow vector S 
is always parallel to the interface and along the posi- 
tive direction of £ axis, implying that the power always 
flows from the gain toward the loss region, and no energy 
exchange between nonlinear and linear media. Because 
1/[1 + e a ^~ r, " ) } 2a is similar to the step function £/(£, 77), 
the power-flow density is mainly concentrated in the non- 
linear medium. 
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FIG. 1. (a) The imaginary and the real (b) parts of the po- 
tential R with a = 5. (c) The profile of (1 — U)AG based on 
Eq. ([lip . (d) and (e) The intensity distribution of the analyt- 
ical surface solitons with a = 1 and 5, respectively, (f) The 
maximum intensity of the propagating surface solitons versus 
the propagation distance with 5% noise for different a. The 
dot-dash line displays the interface between the two media. 
The other parameter is Wo = 0.1. 



The imaginary and real parts of the potential i?, and 
(1 — U)AG based on Eq. ([TT]h respectively, are displayed 
in Fig. [[] (a), (b) and (c) with a = 5 and W = 0.1. 
The intensity profiles of the analytical surface solitons 
are plotted in Fig. [1] (d) and (e) with a = 1 and 5, re- 
spectively. It is clearly seen that the center of the soliton 
profile does not locate at the interface when a is bigger, 
but is slightly away from the interface into the nonlinear 
medium along the odd symmetric axis (£ + 77 = 0) of the 
VT symmetric potentials so that only a small fraction 
of the soliton power is existed within the linear medium, 
as shown in Fig. [1] (e), which is the result of the initial 
net transverse power flow from the linear to nonlinear 
media. The intensity profile of the soliton is symmetrical 



in both gain and loss regions. Such surface solitons can 
be formed and stabilized because the net gain and loss, 
respectively in the gain and loss regions is zero due to 
transverse power flow from gain to loss regions, while the 
power flow arises from the nontrivial phase structure of 
the surface solitons [9[, as the soliton solution shows in 
Eq. (p2j) . We have checked the robustness of the analyt- 
ical surface soliton by simulating the soliton propagation 
with 5% random noise imposed on both amplitude and 
phase. The maximum intensity of the soliton in propa- 
gation is plotted in Fig.[T](f), showing the stability of the 
soliton for different a. One can see that the stability of 
the soliton is very well. 




FIG. 2. (color online) The intensity profile of the solitons 
with the propagation constant b — 2.44 (a) and 7 (b). (c) 
The power P of the surface soliton versus b. The red and 
black filled circles correspond to the solitons in (a) and (b), 
respectively, (d) The real part of the perturbation growth 
rate 5 versus b. (e) Transverse power flow of the soliton in 
(a) where the contour plot in the background represents the 
its intensity distribution, (f) The maximum intensity of the 
surface solitons in (a) and (b) propagate with 5% noise. The 
other parameters include Wo = 0.1, AC = — 1. 

From Fig. [1] (d) and (e), one can see that the power 
of the surface solitons confined in the linear medium is 
very small or almost zero when |AG| or a is bigger. Ac- 
cordingly, we can draw a conclusion that VT surface 
solitons can exist and stabilize even if VT symmetric 
potential exists only in nonlinear medium regardless of 
VT symmetry of the linear medium because only very 
small portion of the soliton power is in linear medium. 
To check the conclusion, the VT potential in Eqs.JSJ)- 
(fTTj) is modified as V = [/"(£, r/) [(4 + ^f) sech 2 (£ + rj) + 
4sec/i 2 (£ — rj) — sech 2 (^ + rj) sec h 2 (£ — rj)], and W = 
WoU(£, rj) sec h(£ + rj) tanh(£ + rj), that is to ignore the 
VT potential in the linear medium. Numerical soliton 
solutions to Eq.flBJ) under the modified VT potential are 
obtained using the spectral renormalization method [19| , 
and are displayed in Fig. [2] (a) and (b) with b = 2.44 
and b = 7, respectively, for AG = —1 and Wq = 0.1. 



The spot size of the solitons is obviously affected by the 
propagation constant b, as shown in Fig. EJa) and (b). 
Furthermore, the intensity of the solitons increases with 
b. The power of solitons, P = J_ f_ \u\ 2 d^drj, is plot- 
ted in Fig. [2] (c) as a function of b, and is obviously lower 
than the self-focusing critical power (P cr « 11.70) [21|. 
As a consequence, the formation of the VT surface soli- 
tons is no requirement of critical power threshold. 

The stability of the surface solitons is also checked with 
the linear stability analysis method. It is assumed that 
<ft£, V) = «(£, V)e ib( + e[F(£, V )e S( + G*(£, r?)e^]e*« [9, 
[20[, where e « 1, F and G are the perturbation eigen- 
functions and 5 is the growth rate of the perturbation. 
Substituting q'{^,rj) into Eq.([7]) and linearizing the equa- 
tion, we gain, 

5F = z[(Vi + R + i) + 2U\u\ 2 )F + Uu 2 G], (14) 
5G = -i[(Vi + R*u + i) + 2U\u\ 2 )G + Uu* 2 F], (15) 

where ^(£,77) = [1 - U(£,r])]AG - b. 

It is well known that the soliton is linearly unstable if 
the real part of growth rate S is positive. Otherwise, it is 
linearly stable [20|. Fig. [2] (d) shows the Re(S) versus b. 
One can see the real part of S is not positive as b > 2.35. 
As a result, the surface solitons are always stable until the 
power of the solitons is close to the self-focusing critical 
power P cr . Numerical results also show that the surface 
solitons do not exist as b < 2.35. The transverse power 
flow (Poynting vector) is mapped by arrows in Fig.[2](e), 
which intuitively indicates the energy exchange among 
gain or loss domains. The robustness of the surface soli- 
tons has also been checked with 5% random noise. The 
maximum intensity of the solitons in propagation is plot- 
ted in Fig. [2] (f), showing the stability of the solitons for 
different b. 

Next, setting AG = — 50 (corresponding to a larger 
refractive difference Uq 2 — ^oi m ^ ne modified potential), 
numerical surface solitons are calculated and plotted in 
Fig. [3] (a), (b) and (c) for different propagation constant 
b. One can see that the power of the surface soliton 
confined in the linear medium is obviously less than one 
in Fig. [2] (a) and (b) where |A£?| is smaller than here. The 
total power of the solitons can also be found still lower 
than the self-focusing critical power, as Fig. [3] (d) shows, 
and hence no power threshold required for the formation 
of surface solitons. Fig. [3] (e) shows the stable region of 
the surface solitons, which is similar to one in Fig. El (d). 
The maximum intensity of the solitons in propagation 
is plotted in Fig. [3] (f), also showing that the soliton is 
stable, and revealing the agreement with the analysis of 
linear stability in Fig. [3] (e). 

To explore the generalization of the VT surface soliton, 
its existence and stability are investigated in another kind 
of potentials, a two-dimensional Scarff II potential [9[, 
which is expressed by R = U(£,rj){Vo sec h 2 (^/^ 2 + rj 2 ) + 
iW sech(^£ 2 + 77 2 )[tanh(£) + tanh(^)]}, where Ufarj) 
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FIG. 3. (color online) The intensity profile of the solitons 
with b = 2.44 (a), 3 (b) and 7 (c). The dot-dash line displays 
the interface between the two media, (d) The power P of the 
surface solitons versus b. The red, black, and pink filled circles 
correspond to the solitons in (a), (b) and (c), respectively. 

(e) The real part of the perturbation growth rate 5 versus b. 

(f ) The maximum intensity of the surface solitons propagates 
with 5% noise. The red, black and pink curves correspond 
to the solitons in (a), (b) and (c), respectively. The other 
parameters are Wo = 0.1, AC = — 50. 
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FIG. 4. (color online) (a) The real and (b) the imaginary 
parts of the modified Scarff II potential R. The intensity 
profile of the solitons with b = 0.16 (c), 0.84 (d) and 1.46 (e). 
respectively, (f ) The power P of the surface solitons versus 
propagation constant b. The red, black and pink filled circles 
correspond to the solitons in (c), (d) and (e), respectively. The 
dot-dash line displays the interface between the two media. 
The other parameters are Vq = 3, Wo = 0.3, AC = — 1. 



is still a step function as depicted before. The numerical 
VT surface soliton solution to Eq.flSJ) is obtained under 
the two-dimensional Scarff II potential. The real and 
imaginary parts of the potential R are displayed in Fig. 2] 
(a) and (b), respectively, while the typical surface solitons 
are plotted in Fig. 2] (c), (d) and (e), respectively with 
b = 0.16, 0.84 and 1.46 for V = 3, W = 0.3 and AG = 



— 1. The power of the soliton is plotted in Fig. [4] (f) as 
a function of 6, and is again shown lower than the self- 
focusing critical power P cr , which is similar to the result 
in above other potentials. 

The stability of surface solitons may be measured by 
the real part of perturbation growth rate. The real parts 
of the growth rate versus b are plotted in Fig. 0(a). One 
can find that the stability of the solitons maintains only 
in a small range of 6, which is much different from the 
case in the other potentials before where the soliton can 
stabilize in a wide range of 6, as shown in Fig. [2] (d) 
and Fig. [3] (e). Fig. [5] (a) shows the soliton in Fig. 0] 
(e) (b = 1.46) is unstable. Fig. 0(b), (c) and (d) show 
the linear stability spectra of the solitons , respectively, 
corresponding to the solitons in Fig. S](c), (d) and (e). 
Fig. [5] (d) still shows the unstability of the soliton in 
Fig. 2] (e). The maximum intensity and power of the 
solitons in propagation are displayed in Fig. [5] (e) and 
(f), respectively. They all show that the soliton with 
b = 1.46 is unstable due to the oscillation of the intensity 
and power in propagation till divergence eventually. 
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FIG. 5. (color online) (a) The real part of the perturbation 
growth rate S versus b. The red, black and pink filled circles 
correspond to the solitons in Fig. [H (c), (d) and (e), respec- 
tively. Linear stability spectra of the soliton with b = 0.16 
(b), 0.84 (c) and 1.46 (d). (e) The maximum intensity and 
(f ) the power of the surface solitons in propagation with 5% 
noise corresponding to Fig. |H(c), (d) and (e), respectively. 



In summary, a novel class of two-dimensional surface 
solitons is reported. It exists at the interface between 
the linear and nonlinear complex refractive index media 



with VT symmetry. The soliton equation governing the 
evolution of the surface soliton is derived. Analytical 
and numerical solition solutions are explored over a wide 
range of VT potential structures and parameters. It is 
found that the surface soliton can stabilize over a wide 
range of potential structures and parameters, and even 
can exist and stabilize at the interface where nonlinear 
medium alone meets VT symmetry but linear medium 
does not. One-dimensional VT surface soliton can not 
exist, which is distinct from the case in single VT sym- 
metric potential where one-dimensional VT soliton can 
exist and stabilize within the VT potential. The cause 
that one-dimensional VT surface soliton can not exist is 
that the gain and loss regions are alone distributed in one 
of the linear and nonlinear media, while the interface be- 
tween the two media prevents free power flow from the 
gain to loss domains. The nontrivial phase structure of 
VT surface solitons plays a key role in maintaining stabil- 
ity of the VT surface solitons. It is the nontrivial phase 
structure that leads to the transverse power flow from 
gain to loss domains to keep energy conservation in the 
propagation of the VT surface solitons along the inter- 
face. Furthermore, these surface solitons do not exhibit 
a requirement of power threshold. 
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